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The space and time evolution of initial parton production in ultrarelativistic 
heavy ion collisions is investigated within the framework of perturbative QCD which 
. includes both initial and final state radiations. Uncertainty principle is used to relate 

o _ 

■ the life time of a radiating parton to its virtuality and momentum. The interaction 
time of each hard or semihard parton scattering is also taken into account. For 

^ . central Au + Au collisions at ^/s=200 GeV, most of the partons are found to be 

■ produced within 0.5 fm/c after the total overlap of the two colliding nuclei. The local 



X 



momentum distribution is approximately isotropical at that time. The implication on 
how to treat correctly the the secondary scattering in an ultimate parton cascading 
model is also discussed. 
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I. INTRODUCTION 



In the search for a quark gluon plasma (QGP), a thermally and chemically equi- 
librated system of quarks and gluons freed from the color confinement, heavy nuclei 
are accelerated to collide with each other at ultrarelativistic energies. Studies in 
the last few years [|I|- suggest that heavy ion collisions are dominated by minijet 
production via semihard parton scatterings with pt ~ few GeV/c at energies of the 
Relativistic Heavy Ion Collider (RHIC) of Brookhaven National Laboratory (BNL) 
and the proposed Large Hadron Collider (LHC) of CERN. These minijets may con- 
tribute to most of the transverse energy produced in the heavy ion collisions. If there 
is enough time and if the density is sufficiently high for the partons to rescatter, the 
initially produced minijets will eventually lead to a thermalized QGP. However, in 
order to estimate the initial energy density of the produced partons, one must know 
their formation time. The total formation time must include the interaction time of 
the semihard processes and the formation time of initial and final state radiation. 

Another motivation of our study in this paper is to understand the space-time 
evolution of the initial parton production and the consequences on secondary scat- 
terings. One has to consider these scatterings in order to simulate the thermalization 
process by parton cascading. Although there are already Monte Carlo simulations 
of parton cascading in heavy ion collisions 0,0, the possibility of double counting is 
still a potential problem. One type of double counting may occur if one neglects the 
interaction time of the scatterings. During the interaction time, which depends on 
whether the scatterings is hard, semihard or soft, the participating and therefore the 
produced partons cannot scatter again immediately with the beam partons. Another 
double counting has to do with the initially radiated partons before the hard or semi- 
hard scatterings. In calculating jet cross sections and the number of hard scatterings, 
parton structure functions of a nucleus fa/A{x,Q'^) evaluated at the scale = pip 
are used. These parton distribution functions then have already included the effect 
of QCD evolution, producing more partons at small x and larger . Therefore, in 
Monte Carlo simulations in which initial radiations are treated as backward evolu- 
tion, the radiated partons should not participate in any interactions before the end of 
the corresponding hard scattering. This is also in coincidence with the factorization 
theorem [|| of perturbative QCD (pQCD). 

In this paper, we will use HIJING Monte Carlo model [Q], which was developed 
for simulating parton and particle production in heavy ion collisions, to study the 
space-time evolution of the initial parton production. In HIJING, each hard parton 
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scattering in hadronic reactions is simulated by using some subroutines of PYTHIA 
1^ in which one can trace back the whole history of initial and final state radiation 
associated with that hard scattering. The life time of an intermediate off-shell parton 
is estimated via the uncertainty principle tq^ Ejc^ for given energy E and virtuality 
q. We then obtain the space and time vertices of all the produced partons. Given 
a space cell, we then calculate the time evolution of the momentum distribution of 
the produced partons by assuming classical trajectories without rescatterings. We 
will verify that after some time most of the partons will leave the space cell except 
those whose rapidities are close to the spatial rapidity of the cell, thus achieving local 
isotropy in momentum distribution which happens to look like thermal. This will 
provide a good starting point for treating equilibration via a set of rate equations 
as in Refs. for the number and energy densities. However, what is not studied 
in this paper, though included in HIJING, is the particle production from the soft 
component. 

This paper is organized as follows: In Sec. I, we briefly review the HIJING model 
and the related aspects of parton shower in PYTHIA. More detailed descriptions can 
be found in Refs. In Sec. II, we discuss how we implement the Lorentz contracted 
spatial parton distribution of a nucleus into the calculation and how we estimate the 
life time of a virtual parton and the interaction time of semihard processes. The results 
are presented and discussed in Sec. HI, where consequences for secondary parton 
scattering and thermalization processes are also discussed. We give our conclusions 
in Sec. IV. 



II. HIJING MODEL 



To make a consistent study of the initial parton production, we briefly review the 
HIJING model and the structure of initial and final state radiation in PYTHIA which 
HIJING model utilizes exclusively for hard parton scatterings and the associated 
radiations. Readers can find detailed descriptions in Refs. 



A. Minijet Production 

In perturbative QCD, given the parton structure functions fa/N{x,Q^) and the 
perturbative parton-parton cross section daab, the differential cross section for jet 



production in nucleon-nucleon collisions can be calculated via [jTO| 
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fa/N (Xi ,pI)x2 fb/N ,pI)-^ 



dp^dyidy2 ^ dt 

cd 



where Xi^2 are the fractional momenta of the colhding partons, 7/12 are the rapidities 
of the jets and Pt is their transverse momentum. The factor K ^ 2 accounts roughly 
for the corrections beyond the leading order. We want to emphasize here that the jet 
cross sections are calculated with the parton structure functions evaluated at the scale 
of the hard scattering = p^. These structure functions fa/N{x, Q"^) are obtained by 
evolving them from a lower scale Qq via Altarelli-Parisi equations [|ll|] thus increasing 
the parton density at lower values of x. 

Given the integrated inclusive jet cross section ajet, we can use the unitarized 
eikonal model [ll^l- IITTII to calculate the total nucleon-nucleon cross section and the 



Tjv(b) 



cross sections for multiple independent jet production |18] 

ao = J d% [1 - e-"-"^^(^) 



^2^het^n^g-.,..T,(6)^ (3) 



where T/v(&) is the partonic overlap function between two colliding nucleons at impact 
parameter b. 

In this model two phenomenological parameters have to be introduced |TE[: An 



infrared cutoff Pq = 2 GeV/c is used to calculate the total inclusive "hard" parton 
interaction cross section crjet(po) V^)- Another parameter, cTsoft = 57 mb, is used to 
characterize the corresponding "soft" parton interactions. Most of the differences 
among several existing models 0-0] result from different values of these two param- 



eters used. As demonstrated in Ref. |18[, these two parameters, though constrained 
by the total pp and pp cross sections, are still model dependent. Since most of the 
minijets are nonresovable as distinct hadronic clusters in the calorimeter of an exper- 
imental detector, their existence can only be justified by their contribution to particle 
production in low and intermediate pt region. Therefore, the value of po and cTgoft 
will depend on what is included in the particle production from the so-called "soft" 
interactions. Recently, two-particle correlation functions in azimuthal angle and their 
energy and px dependence are proposed to give further constraints on these two 
parameters. 

Extrapolated to heavy ion collisions which are decomposed into binary nucleon- 
nucleon scatterings, Eqs.(|^J^) are used to determine the number of minijet production 
in each binary collision. Each hard scattering is then simulated via subroutines of 
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PYTHIA and initial and final state radiations are generated. At ultrarelativistic en- 
ergies, the colliding nuclei are highly contracted longitudinally in the center-of-mass 
frame but are surrounded by soft partons with the same longitudinal size as a nucleon. 
The hard or semihard scatterings all happen at around the same time when the two 
nuclei pass through each other. Therefore, binary collision is a good approximation 
for hard scatterings. In each binary collision the same parton structure can be used. 
In other words, each individual nucleon really does not have time to readjust its par- 
ton distribution before the next binary collision. Of course, in the actual simulation, 
energy and momentum conservation is strictly imposed. To take into account the 
nuclear modification of the parton structure functions, we parameterize the nu- 
clear shadowing effect as measured experimentally in deeply inelastic lepton-nucleus 
scatterings and assume gluons and quarks are shadowed by the same amount in small 
X region. The inclusion of nuclear shadowing will effectively reduce the number of 
minijet production by about %50 at RHIC energy. However, considering the QCD 



evolution of the shadowing |20] may reduce the effective shadowing. 



B. Initial and Final State Radiation 

For each hard scatterings, one then has to take into account the corrections due to 
initial and final state radiations. In an axial gauge and to the leading pole approxima- 
tion, the interference terms of the radiation drop out. The amplitude for successive 
radiations has then a simple ladder structure and the probability for multiple emission 



becomes the product of each emission [^. The virtualities of the radiating partons 



are ordered along the tree, decreasing until a final value /Iq below which pQCD is no 
longer valid any more. This provides a framework for a Monte Carlo simulation of 
parton shower and its space-time interpretation p2| , ^ . 

At a given vertex of the branching tree, the probability for the off-shell parton a 
of virtuality < q^g_^ to branch into partons h and c with fractions z and 1 — z oi 
the light-cone momentum is given by p2| , [23[l 



q ZTT ^a\<l ) 



where Pa^bc{z) is the Altarelli-Parisi splitting function |Tl[] for a —> be process. By 



requiring the relative transverse momentum qx of b and c to be real, 

q'T = 4l-z)U-^j-^)>0, (5) 
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and a minimum virtuality g^, > /Iq, the kinematically allowed region of phase space 
is then, 

4/io < dax; 

e(g) < z<l- e(g), e(g) = ^(1 - ^l-Afil/q'). (6) 
The Sudakov form factor iSa(g^) is defined as 

y Jifil Je{k) ^ 2-K J 

SO that Sa{q'^s,yr)/Sa{q^) is the probability for parton a not to have any branching 
between g^^^^ and q^. 

In principle, one can perform the initial state radiation in a similar way. The 
partons inside a nucleon can initiate a space-like branching increasing their virtuality 
from some initial value Qq. A hard scattering can be considered as a probe which can 
only resolve partons with virtuality up to the scale of the hard scattering. Otherwise 
without the scattering, the off-shell partons are only virtual fluctuations inside the 
hadron and they will reassemble back to the initial partons. In PYTHIA, which 
uses backward evolution, a hard scattering is selected first with the known QCD- 
evolved structure function at that scale, and then the initial branching processes are 
reconstructed down to the initial scale Qq. The evolution equations are essentially 
the same as in final state radiation except that one has to convolute with the parton 
structure functions. Readers can find details in Ref. 0. 



III. FORMATION AND INTERACTION TIME 

When a parton is off-shell, it can be considered as a virtual fluctuation and it 
can only live for a finite time. At, determined by its virtuality q^ via the uncertainty 
principle. 

At ^ go/g', (8) 

where go is the energy of the parton. After At, the off-shell parton will then branch 
or "decay" into other partons which can further initiate branchings until a minimum 
virtuality fiQ is reached. At g^ < /ig, pQCD is not considered to be valid anymore 
and the process of nonperturbative hadronization takes over. Following this tree of 
branching (which also includes initial space-like radiation) and assuming a straight 
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line trajectory for partons, we can then calculate the space and time evolution of 
the initial parton production. Since we don't consider secondary parton-parton in- 
teractions, the produced final partons which are on shell are considered free particles. 
Fig. |1] is a demonstration of a hard scattering and the associated branching trees. 
The number on each line represents the life time (in fm/c) of that virtual parton and 
the numbers in parentheses are the corresponding virtualities (in GeV). 

The life time of a virtual parton in Eq. |^ is only a rough estimate according to 
the uncertainty principle. One could also estimate it via pQCD in lowest order |Q 
which gives an extra factor of l/a^. But higher order corrections could easily cancel 
this factor, so that q^/q^ should be a good estimate of At in magnitude. 

We also use Eq. ^ to estimate the interaction time for each hard scattering with 
q being the sum of the initial or final four-momentum of the colliding partons. If the 
fractional momenta of the partons are Xi and X2, then the interaction time can be 
estimated as 

2XiX2y/s 

In this case, the asymmetric scatterings (xi ^ X2 or xi ^ X2) have longer interaction 
time than the symmetric ones (xi ~ X2) for fixed parton-parton center-of-mass energy 
Xi X2\fs- We also assume that the interaction time is the same for all channels. 

In the rest frame of each nucleus, three-parameter Wood-Saxon nuclear densities 
are used to construct the nucleon distribution inside the nucleus. The system is then 
boosted to the center-of-mass frame of the two colliding nuclei. Due to the fact that 
gluons, sea quarks and antiquarks are only quantum fluctuations before they really 
suffer scatterings, their longitudinal distribution around the center of the nucleon is 
still governed by the uncertainty principle in any boosted frame [^. We refer to this 
distribution as the "contracted distribution", in which a parton with Xi fractional 
momentum has a finite spatial spread, 

A;zi ^ Ijxi^s. (10) 

Transversely, partons are distributed around their parent nucleons according to the 



Fourier transform of a dipole form factor |18|. If we define t = as the moment when 



the two nuclei have complete overlap, then the interaction point of two partons in a 
t — z plane can be anywhere within the shaded area in Fig. |[ The solid lines are the 
trajectories of the two parent nucleons which spread around the nuclei according to 
a longitudinally contracted Wood-Saxon distribution. 
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Because of the spread of parton interaction points in space and time, it is am- 
biguous to define a common proper time in the t — z plane for parton production in 
heavy ion colhsions. Therefore, as we will demonstrate in the next section, a proper 
time, 

T=Vt^- Z\ (11) 

defined with respect to the nucleus overlap point, is only approximately relevant. 



IV. NUMERICAL RESULTS 

In this section, we present the numerical results of our study. All our results are for 
central Au + Au collisions at RHIC energy. In our simulation, Duke-Owens parton 



distributions p5[ for nucleons are used. The initial virtuality for the initial state 
evolution is set to be Qq =2 GeV/c, and the minimum virtuality for the final state 
radiation is fiQ = 0.5 GeV/c. The maximum virtuality for the associated radiations in 
a hard scattering with transverse momentum transfer is chosen to be gmax = '^Pt- 
In PYTHIA, angular ordering is also enforced to take into account the soft gluon 



interference p3[ in the final state radiation. In HIJING model, soft interactions as 



string excitations are also included. These soft interactions must also consume energy 
and affect minijet production slightly. We want to emphasize that in our calculation 
here, we do noi include the soft interactions, so the number of minijet production here 
should be slightly more than the full calculation when soft interactions are included. 



A. Space-time Evolution 

To estimate the parton production time, we plot in Fig. ^ the total number of 
produced partons, on-shell as well as off-shell, as a function of time at the highest 
RHIC energy, ^/s = 200 GeV/n. We see that long before the two nuclei overlap 
and hard scatterings take place, partons have already been produced via initial state 
bremsstrahlung. Some of the initially radiated partons will also initiate time-like 
branching trees(see Fig.|l|). Note that, if the coherence is not broken by the hard 
scattering, partons which would have been emitted from the initial state radiation 
will not emerge as produced partons. Here we have also included the initiators of 
the space-like branching as produced partons. Therefore, if a parton does not have 
initial state radiation, it will only become a produced parton after the hard scattering. 



8 



whereas, a parton is defined to be produced before tlie liard scattering if it has initial 
state radiation. From Fig. ^, we can see that about 2/3 of the total number of partons 
are produced between t = —0.5 and 0.5 fm/c while about 200 semihard scatterings 
happen between t = —0.1 and 0.1 fm/c as indicated by dashed lines. We find also 
that about 2/3 of the total number of partons are produced in initial and final state 
radiations. The fraction of partons from branching should increase with the colliding 
energy and with smaller choices for fiQ. 

To see how hard or semihard scatterings and initial and final state bremsstrahlung 
contribute to parton production, we show in Fig. § the rapidity distribution of pro- 
duced partons at different time. Before t = 0, most of the partons come from initial 
state radiation. Since the radiations are almost collinear, these partons move along 
the beam direction and therefore have large rapidities. The semihard scatterings then 
produce partons uniformly over a rapidity plateau and fill up the middle rapidity re- 
gion. Final state radiations, which happen after the semihard scatterings, will also 
produce partons uniformly in the central rapidity region. We therefore see from Fig. ^ 
that the dip of dN/dy in middle rapidity is actually caused by the parton produc- 
tion from initial state bremsstrahlung at large rapidity. These partons with large 
longitudinal momenta will move away from the interaction region after the semihard 
scatterings. They do not rescatter with the beam partons in the leading twist approx- 
imation. We will come back to this point again when we discuss the consequences 
of our study on how to treat parton rescattering. In Fig. ^, we also show the time 
evolution of the pr distribution of the produced partons. Since the partons, which 
initially have a Gaussian distribution, have gone through initial state radiation, 
they have already a large pt tail at t = —0.6 fm/c before the hard scattering. Hard 
scatterings will transfer large transverse momentum to the final partons and the time- 
like branchings produce a lot of partons with small Pt- The final Pt spectrum then 
looks more or less like an exponential one, even though secondary scatterings have 
not yet been taken into account. 

In Fig. we plot dN / dz as a function of z at different time to illustrate how the 
parton production evolves in space and time. At t = —0.7 fm/c, as the two nuclei 
approach toward each other before they actually overlap, initial state radiations have 
already begun. These partons have large rapidities and are Lorentz contracted with 
an average spread in z, 

Az ^ 1/po + 2i?A^^ ~ 0.25fm, = 200GeV, (12) 

where Pq = 2 GeV/c is the px cutoff for semihard scattering, Ra is the nuclear radius 
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of Au and is the nucleoli mass. After the hard scatterings and during the interac- 
tion time, partons are produced uniformly in the central rapidity region. Afterwards, 
partons follow straight line by free-streaming and are distributed evenly in z between 
two receding pancakes of beam partons (partons from initial state radiations). One 
can clearly see that there is approximately boost invariance in the central region in 
this free-streaming picture. 

Another illustrative way to study the evolution of local parton density is to make 
a contour plot of parton density p in z and t as shown in Fig. 0. Here, p is defined as 

1 dN 
7tR% dz 

and a sharp sphere distribution is assumed for nuclear density. One can clearly see 
that partons inside the two approaching nuclei have a spatial spread of /S.z = 0.25 
fm in z. This spread continues for partons from the initial state radiation as they 
escape from the interaction region along the beam direction with large rapidities. The 
interaction region where semihard scatterings happen lasts for about 0.5 fm/c, from 
t = —0.25 to 0.25 fm/c. Because of the contracted distribution, some of the partons 
lie outside the light-cone which is defined with respect to the overlapping point of the 
two nuclei. The definition of a proper time r (Eq. |ll]) for the evolution of the whole 
system is therefore only valid within an accuracy of Az. 



If one assumes boost invariance ||26|, the parton density can be estimated as 

irKj^T dy 

where r is the proper time and p should be a function of r only. By comparing the 
contour of constant density with the hyperbola of constant r (dot-dashed line) in 
Fig. 1^, we see that this is true only approximately. We also see that the density, as 
indicated by the numbers, decreases like 1/r due to free-streaming. 



B. Local Isotropy in Momentum Space 



Since there are numerous partons produced within a rather short time as we have 
demonstrated, the initial parton density is very high at t > 0.25 fm/c, immediately 
after the interaction region (see Fig. |^). Within such a dense system, secondary parton 
scatterings and production are inevitable. The equilibration time for the system can 
be estimated by solving a set of rate equations as recently has been done in Ref. [||. 
In this approach, one must make sure that there is approximately local isotropy in 
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momentum space. This could be achieved through secondary parton scatterings as 

has been investigated in Refs. PJ27[]. One can also use free-streaming to estimate 

the upper bound of the thermalization time tiso P,p8|, by studying the momentum 

distribution of partons in a cell of the size of the mean-free-path A/. We take A/ ~ 1 

fm in our study here. 

Let's concentrate on the central slice at z = with \z\ < 0.5 fm. At the very early 

stage during the interaction region (for example at t = in Fig. |^ ), produced partons 

with different rapidities are confined to a highly compressed slab with Az ^0.5 fm. 

As the system expands, partons with large rapidities will escape from the central 

slice while partons with small rapidities remain. As shown in Fig. ||, the rapidity 

distribution of the produced partons in this central slice evolves from a plateau-like 

distribution with width Ay/ 2 ^ 2.5 at early times to a Gaussian shape at t = 0.7 

fm/c. If the free-streaming continues, the rapidity distribution in this central slice 

will eventually become a (5-function and only partons with zero rapidity remain. The 

evolution in other cells is similar in their local frames which have spatial rapidity 

1 ^ t + z , , 

r, = -ln^ (15) 

In Fig. 1^, we show the evolution of momentum distributions in (solid lines) 
and Pz (dashed lines) at different times (indicated by the number on each line). The 
Px distribution should evolve like pt distribution as is shown in Fig. ^ The slope of 
Pz distribution, however, decreases because partons with large longitudinal momenta 
gradually escape from the central region. At t = 0.7 fm/c, the slopes of p^ and Pz 
distributions become the same. If one wants to use the parton production we have 
calculated so far as an initial condition to study thermal and chemical equilibration, 
this should be the starting point where local isotropy is approximately achieved. 
Therefore, the time for achieving isotropy in momentum space via free-streaming is 

tiso^O.7 fm/c. (16) 

This is the same as estimated in Ref. |Q . We emphasize that this is only an estimate of 
the upper bound. However, without secondary scatterings, thermalization can never 
be achieved and maintained. 



C. Chemical Composition 

Unlike an ideal gas of quarks and gluons in which chemical equilibrium is main- 
tained, the initial production of quarks and gluons is determined by the parton struc- 
ture functions, the hard scattering cross sections and the radiation processes in pQCD. 
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Due to the difference in the numbers of degrees of freedom in color space, cross sections 
involving gluons are always larger than those of quarks. For small angle scatterings, 
one can verify that, 

^^a^, (17) 
dt ^ ^ 

where t is one of the Mandelstam variables and, 

4 9 

a=g,l,^, (18) 

for i = qq, gq and gg scatterings. Similarly, both initial and final state radiations 
produce more gluons than quarks and anti-quarks. Therefore in pQCD, the ratio 
between produced quarks and gluons is much smaller than the ratio of an ideal gas. 



which is 9/4 for three quark flavors. Shown in Fig. |T0| are the fractions of produced 
quarks and anti-quarks as functions of time t in pp and AA collisions at -^i = 200 
GeV/n. For Au + Au collisions, hard scatterings with pr > 2 GeV/c produce about 
13% quarks and anti-quarks. If initial and final state radiations are not included 
the ratio jumps to 18%, because radiations produce more gluons than quarks and 
antiquarks. In pp collisions, about 28% of the partons produced via hard scatterings 
without radiations are quarks and anti-quarks. The difference between Au + Au 
and pp collisions is due to the different A dependence of the valence and sea quark 
production. Since the partons are produced via binary collisions, the number of 
produced gluons and sea quarks and antiquarks scales like A'^^^. On the other hand, 
baryon number conservation requires valence quark production to scale like A. Given 
the fraction of total quark qpp = 0.28 and valence quark Vpp = 0.14 production in pp 
collisions, one can find for Au + Au collisions, 

Qpp — (1 — ^ ^^^)Vpp /-.qN 

"'^ = 1 - (1 - A-V3).„ ■ 
which gives qAA = 0.18 for A = 197 as we obtained from Fig. ^ of the numerical 



calculation. To demonstrate the effective valence quark production, we plot in Fig. |TI 



the rapidity dependence of the fractional quark number. In pp collisions, the valence 
quark production peaks at large rapidities. In nucleus-nucleus collisions, the valence 
quark production has a different scaling in A than gluons and sea quarks. This is why 
the fraction of quark and antiquark production is suppressed more at large rapidity 
than in the central rapidity region. The relative quark and antiquark production in 
the central rapidity region is further suppressed by final state radiation. 

The small fractional quark production within pQCD has important consequences 
for chemical equilibration of the partonic system. Because of the small initial relative 
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quark density and small quark production cross section as compared to gluon produc- 
tion, it takes a very long time, if ever, for the system to achieve chemical equilibrium. 
If this time is longer than the phase transition time, a fully equilibrated QGP may 
never be formed [Q. 



D. Consequences on Parton Rescatterings 

So far we have not considered final state interactions among the produced partons. 
Though an exact quantum field treatment of parton cascading is not possible with 
present technology, various semi-classical approximations has been undertaken 
The most important and difficult effort in these semiclassical approximations is to 
emulate the quantum effects, like interference, coherence, and especially, nonpertur- 
bative phenomena in QCD. Since we have studied the initial parton production, a 
discussion of the consequences on how to treat final state parton rescatterings seems 
to be in place now. 

We have explicitly taken into account the interaction time for the semi-hard scat- 
terings, which is roughly tj ~ ^/vt- Inside a highly Lorentz contracted nucleus, the 
spatial spread for partons which could participate in a hard scattering with transverse 
momentum transfer pt is also about Az ~ 1/pt- This would leave the produced par- 
tons no time to have another hard scattering of px with the incoming beam partons. 
For finitely contracted nuclei at relatively low energies, this kind of double high pt 
scattering is still possible. However, these kind of higher twist processes should be 
suppressed by a factor oi l/p\. It is also possible for a parton to go through a hard 
and a soft scattering subsequently since the soft partons always have a spatial spread 
of 1 fm. This kind of hard-soft multiple interactions constitute the leading contribu- 
tion to higher twist corrections to hard processes in nuclear collisions p9|. However, 



double semihard scatterings at high energy with px ~ po will be suppressed due to 
finite interaction time. 

The large pt enhancement of both Drell-Yan dilepton production and single 
hadrons in pA collisions at ^/s < 50 GeV is considered as a result of multiple parton 
interactions At these energies, the contracted length of a heavy nucleus is still 
relatively larger than the interaction and formation time. The factorized form of par- 
ton model is then modified due to the finite beam energy. However, a collection of 



experimental data [0,0 even in this intermediate energy range has already shown 
the effects of the finite interaction time. As the energy increases, the interaction 
time becomes more important as compared to the size of an increasingly contracted 
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nucleus. The partons then have less time for secondary scatterings. This then leads 
to the observed decrease of the large pt enhacement in pA interactions as the beam 
energy increases. At ultrarelativistic energies, one should therefore expect the sup- 
pression of double semihard scatterings. 

In this study, we used a QCD-based probabilistic model to calculate the number of 
multiple independent semi-hard parton scatterings. When calculating the inclusive jet 
cross section in Eq. we used the QCD-evolved parton structure function fa/A^x, Q^) 
at the scale of the hard scattering Q^. Therefore by sampling the partons from a joint 
distribution function of two colliding nuclei at Q^, 

fa/A{xi,Q')h/B{x2,Q'), (20) 

we have already included the probability for partons from initial state radiation to par- 



ticipate in the semihard scatterings. As illustrated in Fig. |T^, the two partons which 
are chosen to have hard scatterings may originate from the same initial branching tree. 
Since not every parton has to have hard scattering in our probabilistic model, other 
partons produced from the initial state radiation tree have been chosen not to suffer 
hard scatterings. These partons should not be allowed to have further scatterings 
with the beam partons, even though they are produced long before the hard scatter- 
ing and have to "go through" the interaction region during the space-time evolution 
(see Fig. |^. Another way to understand this is to consider all partons in the initial 
state radiation tree as fluctuations. During the hard scatterings, the beam partons 
act as probes which can resolve the fluctuations up to scale Q^. These partons on 
the branching tree leading to the hard scatterings will become incoherent, while those 
not leading to any hard scattering will eventually reassemble back to the initiating 
partons. In our probabilistic model, the beam partons will be chosen to interact with 
only some of these resolved partons, while others will simply move along the beam 
direction even though they become incoherent to the nucleons. 

In our simulation, partons in different hard scatterings have their own independent 
initial state radiation trees, decreasing the virtuality from Q backward to an initial 
scale Qq below which everything is considered nonperturbative. This corresponds 



to the situation of Q2 = Qo in Fig. |T2|, where both parton fission and fusion are 
included in the nonperturbative regime. However, the two branching trees could also 
be correlated as in Fig. |l^if Qi = Qq. In this case, parton fusion and fission take place 
in the perturbative regime and the parton structure functions will be modified [p3| . 
One, therefore, should use multiple parton correlation function /(xi, . . . In 
this paper, we simply used a parametrization B to take into account the correlation 
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effect. If this is not enough, we might have overestimated the parton production from 
the initial state radiation. 



V. CONCLUSIONS 

We have studied the space-time evolution of initial parton production in ultrarel- 
ativistic heavy ion collisions. At RHIC energy, we found that the production time 
is about 0.5 fm/c after the two colliding nuclei have complete overlap. At t = 0.7 
fm/c, the produced partons inside a cell \z\ < 0.5 fm have momentarily achieved local 
isotropy in momentum space by free-streaming, consistent with the estimate in Ref. 

This could provide a starting point if one wants to investigate the equilibration of 
the produced partons via rate equations. The initially produced partons are however 
far away from chemical balance due to small quark production cross sections. 

The formation time of the produced partons has serious consequences on the 
subsequent secondary scatterings. The number of double semihard scatterings should 
be suppressed if one takes into account the formation time. Due to factorization, 
partons from the associated initial state radiation will also not scatter again with the 
incoming beam partons. 

In hadronic collisions, minijets become important only at above ^/s = 50 GeV/c. 
Soft interactions, however, still have important contributions to the total cross section 
and particle production even at the highest collider energy presently available. What 
we have calculated in this paper is only the parton production through the semihard 
processes. We have not included the nonperturbative soft processes. In HIJING 
model, the soft interactions have been modeled as string formation, which carries 
energy by means of a color field. Naively, one would think these soft interactions are 
still present in the final state cascading. However, as has been shown in Refs. 0,0, 



the high partonic density produced by the semihard scatterings will screen the color 
field and make the soft component less important. What are left over may only be 
those semihard rescatterings among the produced partons, with the infrared cutoff 
replaced by the screening mass. Then we would use pQCD to simulate the parton 
cascading consistently. We hope the study in this paper will pave our way ultimately 
to such a cascading model. 
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FIGURES 

FIG. 1. Illustration of a hard scattering and the associated branching trees. The values 
of xi and X2 are the fractional energy carried by the in-coming partons. The number for 
each intermediate line is the life-time of that virtual partons At = qo/q"^ (f™); and the 
numbers in the parentheses are the corresponding virtualities q (GeV). 

FIG. 2. The overlap region in space and time of two incoming partons each with spatial 
spread of Azi and Az2, respectively. The interaction point is chosen randomly inside the 
shaded region. Solid lines show the trajectories of the parent nucleons. 

FIG. 3. The total number of produced partons as a function of time t, with t = 
defined as when the two nuclei have complete overlap. 

FIG. 4. The rapidity distribution dN/dy of produced partons at different times t (as 
indicated by the number for each line). 



FIG. 5. pt distributions of produced partons at different times. 



FIG. 6. Parton distribution along the z-axis at different times. 



FIG. 7. Contour plot \n z — t plane of the parton density p of Eq. 13, as indicated by 
the numbers. The wavy structure along the light-cone is only an artifact of the plotting 
program. 

FIG. 8. Rapidity distributions of partons in the central slice of \z\ < 0.5 fm at different 
times. 

FIG. 9. Px and pz distributions at different times for partons in the central slice \z\ < 0.5 

fm. 
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FIG. 10. The fractional number of produced quarks and anti-quarks as a function of 
time for Au + Au (solid) collisions, Au + Au without radiations (dashed), andp + p without 
radiation (dot-dashed) . 

FIG. 11. The final fractional number of produced quarks and anti-quarks as functions 
of rapidity for Au + Au (solid) collisions, Au + Au without radiations (dashed), and p + p 
without radiation (dot-dashed). 

FIG. 12. An illustration of the possible correlation between two initial branching trees 
leading to two hard scatterings. The two branching trees can result from parton splitting 
(fission) and recombination (fusion). The boxes represent hard processes with momentum 
scale Q. 
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